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ABSTRACT 

HD 172189 is a spectroscopic eclipsing binary system with a rapidly-rotating pulsating S Scuti component. It is also a member of the 
open cluster IC 4756. These combined characteristics make it an excellent laboratory for asteroseismic studies. To date, HD 172189 
has been analysed in detail photometrically but not spectroscopically. For this reason we have compiled a set of spectroscopic data to 
determine the absolute and atmospheric parameters of the components. We determined the radial velocities (RV) of both components 
using four different techniques. We disentangled the binary spectra using KOREL, and performed the first abundance analysis on both 
disentangled spectra. By combining the spectroscopic results and the photometric data, we obtained the component masses, 1.8 and 
1.7 M G , and radii, 4.0 and 2.4 R Q , for inclination ( = 73.2°, eccentricity e = 0.28, and orbital period n = 5.70198 days. Effective 
temperatures of 7600 K and 8100 K were also determined. The measured vsin; are 78 and 74 km s -1 , respectively, giving rotational 
periods of 2.50 and 1.55 days for the components. The abundance analysis shows [Fe/H] = -0.28 for the primary (pulsating) star, 
consistent with observations of IC 4756. We also present an assessment of the different analysis techniques used to obtain the RVs 
and the global parameters. 

Key words. (Stars:) binaries: spectroscopic - Stars: fundamental parameters (classification, colours, luminosities, masses, radii, 
temperatures, etc. - Stars: oscillations (including pulsations) (Stars: variables:) 8 Set - Stars: abundances (Galaxy:) open clusters and 
associations: individual: IC 4756 



1. Introduction 

HD 172189 (=BD +5 3864, V = 8.85 mag, a = 18 h 38 m 37.6 s , 
6 = +05 d 27 m 55.3 s , J2000) has the combined characteristics of 
being an eclipsing and spectroscopic binary, pulsating star, and 
member of a cluster (Martin 2003, Martin-Ruiz et al. 2005 — 
MR05 hereafter, Costa et al. 2007, Ibanoglu et al. 2009 — 109 
hereafter). Each of these provide unique constraints that allow 
us to test stellar evolution theories in an independent form: a) 
an eclipsing spectroscopic binary system is fundamental for de- 
termining the absolute global parameters of both stars and the 
system with precision; b) a pulsating star allows us to use the os- 
cillation frequencies to probe the interior of the star, thereby also 
determining the evolutionary state; c) cluster membership has 



the distinct advantage that the properties such as age, metallic - 
ity, and distance can be well-determined. Given the constraints 
imposed by the cluster membership and the binary system on 
the mass, age, and metallicity of the pulsating star, the observed 
seismic frequencies can be used to test and improve the cur- 
rent asteroseismic models. For example, as both components are 
rapidly rotating with periods of 2.50 and 1.55 days, (see Sectj5]i, 
we can investigate the effects of rotation, such as the mixing of 
elements and transport of angular momentum. Several theories 
exist regarding rapid rotation, but as yet, observations have not 
been able to confirm any of these hypotheses. Some examples 
of these unproved theories include understanding the interplay 
between rapid rotation and convective cores, enabling the trans- 
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port of angular momentum both poleward and in the radial di- 
rection (Featherstone et al. 2007), or the existence of an over- 
shoot boundary layer between the convective core and the ra- 
diative region, where mixing of nuclear elements can influence 
main sequence lifetimes (Brun et al. 2004). Such theories can 
be confirmed from detailed seismic modeling, once the global 
parameters of the pulsating star have been determined. 

Since HD 172189 was discovered to be a binary system by 
Martin (2003), several groups have shown a keen interest in this 
object. Dedicated photometric campaigns (Amado et al. 2006, 
MR05, Costa et al. 2007, 109) have begun to reveal the true na- 
ture of this system, and several oscillation frequencies have been 
documented from the time series of the 5 Scuti star. This system 
is moreover a selected target of the asteroseismic core program 
of the CoRoT satellite mission (Baglin et al. 2006a,b, Michel 
et al. 2008), and has been continuously observed from space in 
white light for about 150 days from April to September 2008, 
with the aim of interpreting the pulsations. 

With the prospects of using the observed oscillation frequen- 
cies to study the internal structure of the star, we have compiled 
spectroscopic data taken in 2005 and 2007 from various sources 
with the aim of determining some spectroscopic properties of the 
system, to facilitate the future analysis of this star. We determine 
the radial velocities (RVs) of the individual components using 
various techniques (Sect. [3} to subsequently solve for the orbital 
parameters of the sytem, while also providing an assessment of 
the methods employed (Sect. 0J. We combine the RV data with 
photometric data and present a full orbital and component so- 
lution for this object (Sect. [5]l- We subsequently disentangle the 
spectra (Sect. [6j, estimate the effective temperatures T e g using 
the disentangled spectra, and using synthetic spectra (Sect. |7J 
and perform the first abundance analysis of this object (Sect. [8}. 
Discussion and conclusions then follow. Before beginning our 
analysis, we briefly review the literature of both HD 172189 and 
IC 4756. 

One of the first references to HD 172189 and IC 4756 can be 
found in Graff (1923), where HD 172189 is named star 83 and 
has V - 8.69. Later, Kopff (1943) published an analysis (star 
93, V = 8.86 mag), using a referenced tentative spectral typ- 
ing of A6 from Wachmann (1939). Photoelectric observations of 
IC 4756 were then carried out in 1964 in Lowell Observatory 
(Alcaino 1965), with the purpose of determining the distance 
and the absorption of the cluster as well as to establish a crite- 
rion for membership. This author determined from the measured 
V = 8.73 mag and the colour-colour diagram, that HD 172189 
(Alcaino star 10) most likely was not a member of the cluster, 
however, stated that proper motions would be needed to con- 
firm this. They determined a distance modulus of 8.2 mag corre- 
sponding to 437 pc, and an age of 820 Myr for IC 4756. Herzog 
et al. (1975) subsequently determined the proper motions of 464 
stars in the field of this cluster and estimated a probability of 
89% of HD 172189 (Herzog 205) being a member of the clus- 
ter, while Missana & Missana (1995) obtained a 91% probability 
based on proper motions and the position of the star. 

Schmidt and Forbes (1984) measured a vsin/ of 69 kms -1 , 
where ; is the inclination of the rotation axis (assumed to be the 
same for both stars and equal to the inclination of the orbital 
plane). The spectral type has been somewhat discordant in the 
literature, probably due to its binary nature later discovered by 
Martin (2003). Adding the fact that spectral typing of A stars 
can be difficult and that both components of this (at least) dou- 
ble system (see Sect. |9j are rapidly rotating, it is not a surprise 
that different authors have arrived at several inconsistent results: 
A6 V (Wachmann 1939, Herzog et al. 1975), A7 III (Schmidt & 



Table 1. Summary of spectroscopic observations 



Observatory Instrument Dates #Spectra (Used) 

La Silla 2.2m+FEROS June- July, 2005 17 

OHP 1.52m+Aurelie June, 2005 14(12) 

Catania 0.91m+FRESCO May-Aug, 2005 21 (19) 

CalarAlto 2.2m+FOCES May-June, 2007 11 

La Palma NOT+FIES July, 2007 14(11) 



Forbes 1984), A6 III or A4 III (Dzervitis 1987), A2 V (Costa et 
al. 2007), and A6 (109). Schmidt (1978) measured Stromgren 
photometry of the system: ft = 2.820, (b — y) = 0.258 mag, 
ml = 0.123 mag, cl = 1.055 mag, and estimated E(b - y) = 
0. 15. Using the de-reddenned quantities and the tables from Cox 
(2000), HD 172189 appears to be of spectral type late B or early 
A. The range of these spectral types clearly imposes little con- 
straint on r e jf , luminosity L* and gravity log g. 

With the available photometric data of the eclipsing binary 
system (MR05, Amado et al. 2006, Costa et al. 2007), extra 
constraints can be imposed on some of the fundamental param- 
eters. Very recently, 109 published combined photometric and 
spectroscopic data with estimates of global parameters. They 
suggested that the system has component masses of 2.06+0.15 
and 1.87+0.14 M . The cluster has an age of roughly 1 Gyr 
(Alcaino 1965, Mermilliod & Mayor 1990). The two compo- 
nents are rapidly rotating in a non-synchronous fashion, and the 
orbit is quite eccentric, indicating that the stars are not interact- 
ing and are most likely detached MS stars. Furthermore, the sug- 
gested MS turn-off mass for IC 4756 is 1.8-1.9 M (Mermillod 
& Mayor 1990). 

2. Spectroscopic Observations 

The spectroscopic observations used in this work are sum- 
marised in Table Q] and the following subsections describe the 
observation, reduction, and calibration of the spectra taken at 
the various observatories. All of the spectra were subsequently 
barycentric corrected. 

2.1. Aurelie data 

A total of 14 spectra of HD 172189 were obtained in a time span 
of 10 nights from 15 to 25 June 2005 with the Aurelie spec- 
trograph, mounted on the 1.52m telescope, at Observatoire de 
Haute Provence (OHP), France. The instrument has a grating of 
1800 lines/mm, providing a spectral resolution R 25,000. The 
spectral range used was 4528-4675A/4468^1542A, with expo- 
sure times of 1200/1500, or 3600 seconds. As there is no pipe- 
line reduction available, we used standard IRAF (Tody 1986) 
reduction procedures. The continuum normalisation was per- 
formed manually by fitting a cubic spline. Typical signal-to- 
noise ratio (SNR) values of the spectra are 55-63. 

2.2. FIESdata 

The FIES data were obtained during an observation run of a 
separate object in July 2007 using the 2.5m NOT telescope 
at the Observatorio del Roque de los Muchachos. The Fibre- 
fed Echelle Spectrograph (FIES) is a cross-dispersed high- 
resolution echelle spectrograph. We used the medium-resolution 
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setup with R =* 45, 000. The spectral range is 3640-7455 A, with 
a maximum efficiency of 9% at 6000 A. Separate wavelength 
calibrator exposures of Thorium-Argon (ThAr) were obtained. 
During each of the 5 nights of observations, either 2 or 3 spectra 
of HD 172189 were taken with exposure times of 600 seconds. 
We used the FIEStool reduction software (Stempels, 20043) to 
calibrate the wavelength and reduce the spectra. This tool is opti- 
mised for FIES data, although the reduction is standard and calls 
IRAF to perform some of the tasks. SNR values around 5720 A 
are -70. 

2.3. FRESCO data 

From May - August 2005 a total of 21 spectra of HD 172189 
were observed with the FRESCO echelle spectrograph, attached 
to the 91cm telescope of the M. G. Fracastoro Mountain Station 
at Catania Astrophysical Observatory (CAO), Sicily, Italy. The 
FRESCO spectra, with a resolution R 21,000, span the spec- 
tral range from 4320 to 6800 A recorded on 19 orders. Standard 
IRAF reduction procedures were used and the continuum nor- 
malisation was performed manually by fitting a cubic spline. 
Typical exposure times were 3600 seconds and SNR values near 
5720 A are -45. 

2.4. FEROS data 

The Fiber-fed, Extended Range, Echelle Spectrograph (FEROS), 
mounted at the 2.2m ESO/MPI telescope at La Silla (ESO), 
Chile, has a resolution of R = 48,000 and covers almost the 
complete range between 3500 and 9200 A on 39 echelle orders. 
A total of 17 spectra were taken in June- July 2005 with exposure 
times of between 750 and 900 seconds. We reduced the spectra 
using an improved version of the standard FEROS pipeline (see 
Rainer 2003). We used an automated continuum normalisation 
procedure developed by M. Bossi (INAF OAB-Merate) to nor- 
malise the spectra. Typical SNR values in the region of 5720 A 
are -130. 



B 



1.01 








1.00 






• - - - 


0.99 


L 






0.98 






Phase = 0.93 J 


0.97 






y \ 






0.96 













-300 -200 -100 100 200 300 

Velocity (km s ) 

Fig. 1. Least-squares deconvolution (LSD) profile (open circles) 
calculated at phase 0.93 (based on T min = 2452914.644 HJD 
from MR05) by summing and weighting all of the lines in the 
spectrum. The solid curve is the double-Gaussian function fitted 
to the LSD profile using just the data points within the region de- 
limited by the dotted vertical lines. We have highlighted the ex- 
tracted RVs from the LSD+GAU and LSD+MM methods with 
vertical continuous and dashed lines, respectively. 

1999). The central positions of the Gaussian functions are 
the RVs of each component. 

- LSD+MM: Calculating the first moments (Aerts et al. 1992) 
from the LSD profiles. 

- IRAF: Using the IRAF task FXCOR with a synthetic non- 
broadened template and then using the DEBLEND function 
to compute the RVs. 

- KOREL: Disentangling the spectra using the programme 
KOREL (Hadrava 1995), and determining the RVs by fitting 
the observed spectra with the superposition of the Doppler- 
shifted disentangled spectra. 



2.5. FOCES data 

Observations at the observatory of Calar Alto in Almeria, Spain, 
during 10 nights in May- June 2007 were carried out both in vis- 
itor and in service mode. The full optical range with a resolution 
of -35,000 was recorderd with the echelle spectrograph FOCES 
on the 2.2m telescope. In total, 1 1 exposures of HD 172189 were 
collected each with exposure time of 1200 seconds. The echelle 
data were reduced using the standard IRAF procedures. Typical 
SNR values in the region of 5720 A are -57. 

3. Determination of radial velocities 

As both components of the system are rapidly rotating, the spec- 
tral lines are broadened, making it difficult to identify individual 
lines Doppler-shifted due to the orbital motion. In fact, most pro- 
files are merged and line profile variations due to pulsations are 
clearly present making complicated and delicate the determina- 
tion of the RVs. For these reasons, we used the following in- 
dependent methods to determine the RVs, described in the next 
subsections: 

- LSD+GAU: Fitting a double-Gaussian function to the Least- 
squares deconvolution (LSD) profiles (Donati et al. 1997, 
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3. 1 . Least-squares deconvolution 

To make optimised use of the common information available in 
all spectral lines, we deconvolved several hundreds of individual 
lines by comparison with synthetic masks, into a single 'cor- 
relation' profile with a significantly increased SNR, using the 
Least-squares deconvolution (LSD) method (Donati et al. 1997, 
1999) (cf. Fig. [TJ. The synthetic line masks contain lines from 
the VALD database (Piskunov et al. 1995, Ryabchikova et al. 
1999, Kupka et al. 1999). We used the values of T eff = 8250 
K and \ogg = 4.0 dex (see Sect. [7]) for the spectral template. 
Varying the effective temperature with AT e s = 250 or 500 K 
and the gravity with Alogg = 0.05 or 0.1 dex did not affect the 
central positions of the LSD profiles significantly. The fitted RV 
values obtained from using different templates was used to es- 
timate the error on the RVs. The LSD profiles were calculated 
by taking into account all elements, apart from He and H, in 
the regions 4380-4814 A and 4960-5550 A. This involved de- 
convolving about 3000, 2600, 2200, 1525, or 250 lines for the 
FEROS, FIES, FRESCO, FOCES, and Aurelie spectra, respec- 
tively, and resulted in profiles with SNR of 1100-3300, 500- 
1700, 400-2000, 400-1200, and 200-700. The LSD line profiles 
did not all have a continuum level at 1.0, so they were subse- 
quently normalised, followed by a rescaling to homogenise the 
tcpib/flSESJliteinfaimJith different instruments (see description in 
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Uytterhoeven et al. 2008). The velocity steps of the LSD profiles 
were 2 km s _1 (FIES, FEROS and FOCES), 2.3 km s _1 (Aurelie), 
and 7.8 km s _1 (FRESCO). We note that systematic instrumental 
RV offsets of the order of 1 .5 km s _1 are most likely present (see 
Uytterhoeven et al. 2008). However, due to the small amount of 
datapoints per dataset, we currently are not able to quantify and 
subsequently correct for these instrumental differences. 

3.2. Double-Gaussian fit to the LSD profile 

To determine the RVs of each component we fit a double- 
Gaussian function to each LSD profile. We initially fixed the 
widths of the Gaussians while fitting the central positions (the 
RVs) and the amplitudes, and subsequently allowed all of the 
Gaussian parameters to be fit. Cut-off values in the velocity axis 
were used for each fit because the broad wings of the Gaussian 
profiles do not accurately match the LSD profiles. Fig. Q] shows 
an example of an LSD profile when both components are near 
maximum separation at orbital phase 0.93 (based on reference 
eclipse minimum of HJD 2452914.644 days, and orbital period 
n = 5.70198 days from MR05). The deformations in the profiles 
are most likely due to pulsations, and this inhibits the accuracy 
of the RV. The dotted vertical lines delimit the region that was 
fitted. The solid curve shows the Gaussian fit, and the dashed 
curve shows the rest of this function for further velocity values. 
The model is fit several times using various cut-off values, and 
using several LSD profiles which are calculated from different 
spectral templates. The RVs are defined as the mean values of 
these fits, with the standard deviations defining the errors. The 
RVs corresponding to this spectrum are denoted by the vertical 
continuous lines. This method worked well when the compo- 
nents in the LSD profiles were sufficiently separated, hence RVs 
at conjunction are not available using this method. 

3.3. First normalised moments 

Another tool to obtain RVs for the primary and secondary com- 
ponents from the LSD profiles is to calculate the first normalised 
moments ((v), e.g. Aerts et al. 1992). At quadrature the velocity 
profiles of the primary and secondary stars are well separated. 
In these orbital phases we determined the integration bound- 
aries, within which to calculate (v), from the individual Gaussian 
profiles for the primary and secondary components described in 
Subsect. 13.21 Close to conjunction, the velocity contributions of 
the primary and secondary stars are blended, which complicates 
the definition of the profile boundaries. Therefore, we assumed 
a fixed width of the component profiles, derived from the spec- 
tra in elongation phase. The moments were calculated by deter- 
mining one of the integration borders, and calculating the other 
border assuming a fixed profile width. This method is sensitive 
to the profiles, including any deformations due to pulsations. In 
Fig.Q~]the RVs corresponding to this spectrum are denoted by the 
vertical dashed lines showing an offset of a few km s _1 from the 
values of LSD+GAU. 

3.4. IRAF FXCOR 

FXCOR cross-correlates the observed spectra with a template 
spectrum in the Fourier domain. In our case, the merged (one- 
dimensional) spectra of HD 172189 were cross-correlated with 
a Kurucz synthetic spectrum computed with the approximate 
physical parameters of the components of the binary, i.e., T e g = 
8000 K and 7500 K, logg = 4.00 and solar abundance. 
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Fig. 2. Cross-correlation profiles (CCF) of the observed spectra 
with a template spectrum at two different orbital phases (solid 
curves). The dashed curve is a double-Gaussian fit at phase 0.76, 
obtained by fitting the upper 50% of the profiles. 
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Fig. 3. Fit of the observed spectra by KOREL disentangling. The 
circles are the observations while the solid curves are the super- 
imposed RV Doppler-shifted disentangled spectra (y-shifted for 
clarity). 



The range of the spectrum used for the cross-correlation was 
between 5000 and 6500 A, avoiding the two regions most af- 
fected by telluric absorption lines. The object spectra were fil- 
tered in the Fourier domain with a bandpass filter according to 
the resolution of the data. Most information in the Fourier spec- 
trum is above a certain wavenumber which was computed to take 
into account the resolution of our data. The Fourier transform 
of the spectrum was then multiplied by a ramp function which 
starts rising at wavenumber 10, reaches 1 at 20, starts falling at 
2500 and reaches again at 4500. The cross-correlation func- 
tions (CCF) were then calculated. Two examples of the CCFs at 
different orbital phases are shown in Fig. [2] Gaussian functions 
were subsequently fitted to the CCFs, using the full CCF profile 
(later referred to as IRAFfull), and next using only the upper 
50% of the CCF (later referred to as IRAF), resulting in much 
better central Gaussian fits. 
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3.5. Disentangling of spectra 

KOREL is a FORTRAN based code that fits time series of ob- 
served spectra of a multiple stellar system in the Fourier wave- 
length domain, to decompose them into mean spectra of each 
component, and simultaneously finds the best estimate of the or- 
bital parameters (Hadrava 1995). It also calculates RVs of all 
components by fitting each exposure as a superposition of the 
disentangled spectra. We use the recent version KOREL08 to 
find RVs with a sub-pixel precision (cf. Hadrava 2009). 

As the spectral lines are broad and hence more likely blended 
with nearby lines, we carefully chose several wavelength regions 
where the orbital motion could be clearly observed. We have 
chosen regions containing different numbers of spectral lines: 
4537-4595 A, 493 1-4995 A, 4975-5038A, 5300-5368A, 5506- 
5577A, and 6340-6357A. We only used the data from FEROS, 
FIES, and FOCES because these provided the highest SNR as 
well as allowing sufficient spectral resolution for the combined 
data set. 

Each spectral region was disentangled independently con- 
verging the following parameters: periastron epoch, eccentric- 
ity e, longitude of periastron £2, semi-amplitude of the radial 
velocity curve of the primary component Ka, and mass ratio 
q = K A /K B (=M b /Ma). The orbital period 11 was held fixed 
to the constant 5.71098 days (MR05), because the spectral data 
can not improve this previous determination. Also, note that the 
systemic velocity y is not obtained from disentangling, because 
KOREL does not use any spectrum template. To determine the 
optimal orbital parameters, we inspected both the residual value 
of the KOREL fit as well as the extracted RV curves. However, 
broad spectral profiles enlarge the range of acceptable orbital 
values considered as "good fits". For each of the spectral regions 
we obtain a set of RV measurements. 

Fig.[3]shows the observed spectra (circles) for the spectral re- 
gion 4546—4585 A at orbital phases 0.98 (eclipse) and 0. 1 (sepa- 
rated). The solid curve is the sum of the individual disentangled 
primary and secondary spectra, but x-shifted by their RVs at the 
appropriate orbital phase and, for maintaining clarity, we also 
vertically shift the phase = 0.1 spectrum. 



4. Orbital parameters based on RVs 

From each of the methods explained in Sect. [3] we obtain a set 
of RV measurements. In order to obtain the orbital parameters of 
the binary system we fit each set separately using the standard 
radial velocity equations. The results are summarised in Table [2] 
where the orbital period is fixed at 5.70198 days (MR05), To is 
the offset in days from the defined primary minimum, and x\ 1S 
the reduced^- 2 value. For the LSD+MM method we have SNR 
values instead of observational error measurements, so the x\ 
value is not comparable with the other methods. The uncertain- 
ties quoted are the standard formal uncertainties calculated from 
the fitted parameters and using the observational errors given 
for each RV data point. The results for KOREL are obtained by 
weight combining all of the RV measurements from the individ- 
ual spectral profiles, and subsequently fitting using the standard 
equations, while the uncertainties reflect the variation in the fit- 
ted orbital parameters from the different spectral regions. These 
are the RV data that were later used for the simultaneous pho- 
tometric and spectroscopic analysis (Sect. |5j. We also analysed 
the published RVs from 109, who had used standard IRAF pro- 
cedures to determine these. The results are given under heading 
109. 



We define the primary component as the more massive star. 
The primary is the star that shows the deeper LSD or CCF pro- 
files as well as the line-profile variations due to pulsations (see 
Sect [Hi. 

Fig. [4] shows the phased radial velocity data from each of the 
methods: LSD+MM (crosses), IRAF (diamonds), LSD+GAU 
(squares), and KOREL (triangles — we have included a 
weighted average y shift of -28.02 km s for clarity). The dotted 
curve shows the radial velocity solution given by the parameters 
using KOREL and shifting it vertically by y. 

All of the methods presented return values of e, Q, To and y 
consistent with each other (Table [2}. There is, however, a vari- 
ation among the fitted values of Ka and Kb- This is due to the 
sensitivity of the various methods used to extract the radial ve- 
locities. Two possible reasons that affect the determination of the 
RVs are (1) the effect of the pulsations on the spectral profiles, 
and (2) line-broadening due to rotation making it difficult to ac- 
curately determine RVs (especially at conjunction). 

If we inspect Fig. Q] we see that the double-Gaussian func- 
tion cannot correctly reproduce the shape of the LSD profile. 
In fact, due to the wiggles or variations in the line, the actual 
center of the Gaussian could be shifted slightly to the right or 
the left. We have highlighted the center positions derived using 
the LSD+GAU and LSD+MM methods with the continuous and 
dashed vertical lines. In this figure it can also be seen that the lat- 
ter method is sensitive to the minima of the deepest variations in 
both the primary and secondary profiles. In this case, it results in 
smaller absolute values of the RVs than the former method. By 
inspecting these profiles by eye, it is very difficult to distinguish 
which is the true orbital motion RV. 

For phases between ~0.35 and 0.7 we begin to see blend- 
ing of the component profiles using both CCFs and LSD meth- 
ods. Hence the resulting RVs between these phases will not be 
as accurate. For example, the deviations of the RVs using the 
LSD+MM method from the fitted RV curve in Fig. [4] are ex- 
plained most likely by the sensitivity of this method to the de- 
formations in the line profiles induced by pulsations. 

In Fig. [2] we show the CCF calculated using IRAF proce- 
dures at two different orbital phases. For phase = 0.76 we have 
also plotted the Gaussian fits. Again, it is difficult to extract an 
accurate RV value because the lines are so broad and due to the 
variations (wiggles) in the CCFs. These values were obtained 
using the upper 50% of the CCF (although we draw the full an- 
alytical profile). 

By inspecting the profiles of all of the available spectra (and 
CCFs or LSDs) it is not possible to distinguish between shifts of 
a few km s , and it is these shifts that cause the differences in 
the fitted values of Ka and Kb of ~6 km s _1 shown in Table [2] 

The last method used to extract the RVs was KOREL. We 
choose this as the best solution because it is robust against broad- 
ened/blended lines, it simultaneously fits the orbital solution 
(comparing all of the spectra at the same time), and it is inde- 
pendent of spectral templates. These are the RVs that are subse- 
quently used for the simultaneous photometric and spectroscopic 
light curve fitting. 

5. Combined analysis of photometric and 
spectroscopic data 

In order to obtain a consistent photometric and RV solution for 
this system, we combine the RV data with an extensive photo- 
metric data set (MR05). We include photometry from some more 
recent (2005 and 2007) observational campaigns with the twin 
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Table 2. Spectroscopic orbital solutions based on fitting the radial velocities obtained using different data fitting methods, explained 
in Sect. [3] 
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Table 3. Component and orbital results from simultaneous pho- 
tometric and velocity light curve fitting with PHOEBE. 
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Fig. 4. Phased radial velocity solution (dotted curve) obtained 
using the parameters from Table [2] under the column heading 
KOREL. The RVs obtained by the individual methods are also 
shown. The dashed lines delimit orbital phase at and 1 . 



Danish 6-channel uvby/3 photometer at the 1.5m and the 0.9m 
telescopes at San Pedro Martfr, and Sierra Nevada Observatories 
(data reduced in the same manner as MR05). We also in- 
clude photometric data from the P7 photometer mounted on 
the Flemish Mercator telescope at the Roque de los Muchachos 
Observatory (see Rufener 1964, 1985 for reduction procedure). 
In total there are more than 4000 photometric measurements 
obtained in each of the four Stromgren filters, spanning more 
than 3700 days, and 70 out-of-eclipse data points in the Geneva 
UBBiB 2 VVi filter system. We used PHOEBE (PHysics Of 
Eclipsing BinariEs) (version 0.31a; Prsa & Zwitter 2005) and 
FOTEL (Hadrava 1990) to model these data, however the lat- 
ter was used mainly to confirm the results. PHOEBE is a user- 
friendly software based on the Wilson-Devinney (WD) code 
(2007 v. of Wilson & Devinney 1971). 

5. 1 . Photometric Indices 

Using out-of-eclipse measurements, we calculated the 
Stromgren indices. We obtain (b - y) = 0.262 mag, ml = 0.134 
mag, and cl = 1.058 mag. The de-reddenned indices were 
obtained using the method described in Philip et al. (1976), these 
are E(b-y) = 0.158 mag, (b-y) Q = 0.104 mag, m Q = 0.184 mag, 
and co = 1.026 mag. If we compare these values to Crawford 
(1979) we obtain spectral type A6V or A7V. These values 
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1.78 (0.24) 
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L(L ) 


52.42 (2.85) 
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r eff (K) 
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90.6 (0.06) 




a (Re) 


20.32 (0.06) 




q = K A IK B 


0.96 (0.001) 




i (") 


73.2 (0.6) 





correspond to a single star with a mass range of between 1 .7 
and 1.9 M and a T eS of between 7400-8100 K. The value of 
{b - y)o during the primary eclipse is 0.105 mag, practically 
the same value as the out-of-eclipse value. This implies that 
the temperatures of both stars are very similar, consistent with 
the estimated effective temperature ratio {Tb/Ta) of 1.05 from 
MR05. 



5.2. PHOEBE 

We started our fitting by varying r ei jB between = 8000 and 8250 
K. The value of r e ffA = 7500 (or 7750) K is fixed, according 
to the effective temperature ratio of 1 .05 (MR05) (and in agree- 
ment with the results from Sect. [7j. We also adopted as starting 
spectroscopic values those given in Table [2] and the results from 
MR05 (e = 0.24, i = 73°, and q = 0.90). We searched for the 
best solution by using iteratively as free parameters: the lumi- 
nosity of the primary component La, separation a, q, e, Q, i, and 
y. The radii, masses, and secondary luminosity are subsequently 
derived from these values. 

We modelled the stars as a detached system assuming a ra- 
diative albedo for both stars equal to 1 and gravity darkening 
values of gA — gs — 1 ■ We also attempted to model the data using 
the semi-detached configuration, however, most of the parame- 
ters did not converge. The model parameters from the best fit 
are given in Table [3] The errors are derived using PHOEBE. The 
velocities v sin i are determined from spectroscopy (see Sect. [8) 
allowing the derivation of the rotational periods P, ot . 
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6. Disentangling of spectra 

The disentangled spectra were obtained by imposing the orbital 
solution from KOREL shown in Table [2] The used version of 
KOREL is limited in resolution to 4000 points. However, in most 
cases it was quicker and more stable to work on 2000 points at 
a time, over wavelength regions of about 100 A. For this we had 
to rebin the spectra at the appropriate resolution and carefully 
choose the endpoints of the regions where there was continuum, 
while also ensuring that overlapping spectral regions are present 
to facilitate the merging of all of the wavelength regions. 

We attempted to disentangle the spectra in the wavelength 
range from 4600 to 5800 A because this region is most sen- 
sitive to abundances. However, we only obtained reliable re- 
sults for ~100A-length spectral regions between 4976-5627 A. 
Unfortunately, this implies that a spectroscopic determination of 
T e g from the Balmer lines is not possible from the disentangled 
spectra: a better orbital phase coverage and higher SNR spectra 
are needed in order to disentangle these regions correctly. 

Spectra from binary systems normalised to 1 contain no in- 
formation about the relative light contributions of each of the 
individual spectra. We used the literature values from 109 (La = 
0.61 and Lg = 0.39 where La + Lb = 1), and the fitted luminosi- 
ties given by our composite photometric and spectroscopic anal- 
ysis (La = 0.69 and Lb = 0.31, Sect. [5} to obtain the fractional 
contributions to the total light. We shifted the output spectra to 
continuum 0, scaled the individual spectra by their correspond- 
ing fractional luminosity, and shifted the spectra again to 1 to 
obtain the disentangled individual spectra. 

7. Spectroscopic determination of T eS 

Once the spectra were disentangled, we proceeded to obtain the 
LSD profiles of the individual disentangled spectra (Sect. 13. j} , 
in order to estimate T e ff and logg. We calculated the LSD pro- 
files using a range of templates spanning T e ff (7250-8500 K) and 
\ogg (3.5^1.5), and identified the best values as those that pro- 
duced the best LSD profiles (smoothest and most symmetric). 
Using the luminosity ratios obtained, r e ff of the primary was 
best fit with a template spectrum of 7500 K with logg = 4.0. 
The secondary was best fit with a hotter temperature of 8000 - 
8250 K and \ogg = 3.5, although logg = 4.0 was also a good 
fit. In general, the results are much less sensitive to the value of 
log g than to changes in T e s, 

As the disentangled renormalised spectra depend on the as- 
sumed luminosity ratio of both components, we performed a few 
tests to investigate how these results depend on this imposed 
value. We scaled each of the disentangled spectra arbitrarily, and 
proceeded to identify the best T e g for each of the "synthetic" 
spectra in the same manner as above (as double-blind tests). The 
results were again quite insensitive to the value of log g, while 
the derived values of T e g were consistent with results of 7500- 
7750 K for the primary, and 8000-8250 K for the secondary. 
This test showed that the identified T e g were not obtained from 
the imposed luminosity ratio of the stars, but rather the shape 
of the spectra. MR05 obtained an effective temperature ratio 
(Tb/Ta) of ~ 1.05 (secondary/primary), in agreement with our 
result, while 109 obtained a hotter temperature for the primary 
component than for the secondary (Tb/Ta ~ 0.96). 

We also directly compared composite rotationally- 
broadened RV- and y-shifted Kurucz synthetic spectra to 
the best SNR observations near maximum elongation phase. We 
used templates spanning temperatures of 7000 - 8500 K and 
logg of 3.5 - 4.5 dex. Figure [5] shows the observations (grey) 
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Fig. 5. Composite rotationally broadened RV- and y-shifted 
spectral templates (black) with observations (grey) for various 
combinations of r e jf (K) for the primary and secondary compo- 
nents. 
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Fig. 6. The disentangled observed primary spectrum (open cir- 
cles), and the theoretical fitted spectrum (solid) obtained using 
the abundances from Table H] 

compared to various templates (black) for the temperature- 
sensitive H-y6 line. As shown in this figure, the best fitted 
temperatures are 7500-7750 and 8000-8250 K for primary and 
secondary star, respectively. (The \ogg values used are 3.5). 

8. Abundance Analysis 

We determined the abundances of both components using the 
disentangled renormalised spectra. We compared the results with 
both values of luminosity fraction, and the final abundances var- 
ied slightly, but within the quoted errors. Our analysis follows 
the methodology presented in Niemczura & Polubek (2006) and 
relies on an efficient spectral synthesis based on a least-squares 
optimization algorithm (Bevington 1969, Takeda 1995). This 
method allows the simultaneous determination of various param- 
eters associated with stellar spectra and consists of minimizing 
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Table 4. The abundances of chemical elements for the primary 
and secondary stars. 
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8.45 (0.15) 


8.78 (0.28) 


8.39 
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8.89 (0.13) 
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7.35 (0.16) 
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7.53 
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1/- 


6.17 (0.38) 
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7.81 (0.21) 


7.91 (0.12) 


7.51 
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6.67 (0.21) 




7.14 
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6.28 (0.18) 


6.92 (0.23) 
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3.89 (0.26) 


3.17 
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4.90 (0.18) 


5.38 (0.21) 


4.90 
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4/5 


5.81 (0.14) 


6.31 (0.13) 


5.64 


26 


19/20 


7.17 (0.21) 


7.84 (0.19) 


7.45 


27 
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5.57 (0.02) 




4.92 


28 


6/9 


6.13 (0.18) 


6.89 (0.22) 


6.23 


39 


2/2 


2.72 (0.20) 


2.78 (0.17) 


2.21 


40 


-/I 




3.20 (0.08) 


2.58 



the deviation between the theoretical flux distribution and the 
observed one. 

The atmospheric models used for the synthethic spec- 
tra determination were computed with the line-blanketed LTE 
ATLAS9 code (Kurucz 1993), which handles line opacity with 
the opacity distribution function. The synthetic spectra were 
computed with the SYNTHE code (Kurucz 1993). Both codes 
were ported to GNU/Linux by Sbordone (20050 The identifica- 
tion of stellar lines and the abundance analysis were performed 
on a line list we constructed based on VALE0data(Kupkaetal. 
2000 and references therein). 

The shape of the synthetic spectrum depends on the stellar 
parameters such as T e ff , log g, v sin i, microturbulence RV and 
relative abundances of the elements. The effective temperature 
and surface gravity were not determined during the iteration pro- 
cess but were considered as fixed input ones (T e ff = 7500/8000 
and logg= 4.0/4.0 for primary/secondary components), while ^ 
= 2 km s~ ! was adopted. All of the other above-mentioned pa- 
rameters can be determined simultaneously because they pro- 
duce detectable and different spectral signatures. The theoreti- 
cal spectrum was fitted to the normalised observed one in small 
regions, until the results converged. The obtained chemical el- 
ement abundances, relevant errors with the adopted solar abun- 
dances (Grevesse et al. 2007) are given in Table [4] The errors 
are calculated as cr Tot = [<x 2 b + <x 2 g + o"f ogg ]^ 2 where <r d b is 
the standard deviation of the individual element abundance, crj cS 
and cr\ ogg denote the variation of abundance for AT e ff + 250 K 
and Alogg + 0.2. For those elements where only one line was 
used to estimate the abundance, the error is (cr 2 + cr} W 2 . 

Apart from the errors in T e ff and logg, several systematic 
errors are likely affecting the derived abundances. For exam- 
ple, fixing the restriction of the disentangling technique to 
the spectral range of 4976-5672 A and the imposed luminos- 
ity ratio all contribute to these sytematic errors. In addition to 
this, the adopted atmospheric models and the atomic data (in par- 
ticular the oscillator strengths) also contribute. We are also not 
taking NLTE effects into account, which are crucial for a num- 
ber of spectral lines of several elements (see e.g. Asplund 2005; 

2 These are available online: wwwuser . oat . ts . astro . it/ atmos/ 

3 http : / /ams . astr o . univie . ac . at ~v ald/ 



Fabbian et al. 2006, 2009). However, we have estimated that for 
the neutral C lines used here, the NLTE corrections are negligi- 
ble (private comm. Fabbian). Moreover, even for other chemical 
elements we can safely assume that the differential NLTE effects 
between the primary and secondary stars are small and hence can 
not explain the differences for some of the elements between the 
two stars. 

Figure [6] shows a region of the primary disentangled spec- 
trum (circles) and the theoretical fitting to it (solid curve) which 
is obtained using the abundances from Table [4] The chemical 
composition of the primary star is slightly sub-solar, judging 
from the Fe content, while the secondary star shows enhanced 
abundances (see Fig. [7]). The derived rotational velocities are 
vsini = 78 + 3 km s~'and 74 + 4 km s~'for primary and sec- 
ondary component respectively, based on a comparison between 
the spectra with broadened theoretical profiles. Coupling these 
values with the results for R Q and i from Table [3] we derive the 
following rotational periods: P ro t,A = 2.50 + 0.10 days for the 
primary star and P rot; B = 1.55 + 0.09 days for the secondary. 

We have also performed a comparison abundance analysis 
based on a detailed line abundance approach using the MOOG3 
code. The full procedure is described in Arellano Ferro et al. 
(2001) and Giridhar & Arellano Ferro (2005). The number of el- 
ements analysed is smaller than those shown in Table [4] but the 
overall Fe abundance and the indication that the secondary com- 
ponent is iron rich is consistent with our previous result. This 
analysis yielded the best results when T e s of 7750 and 8000 K, 
and log g of 3.0 and 3.5 were adopted for primary and secondary, 
respectively. 

9. Discussion 

9. 1 . Orbital and System parameters 

From the spectroscopic and the combined photometric and spec- 
troscopic analysis we obtain system parameters that are in agree- 
ment: e = 0.29, Q. ~ 80° and Tq ^ -0.1 (this result is compatible 
with To - 0.0 taking the error into account, and this value indi- 
cates that the resulting offset in days from eclipse/conjunction 
coincides with the input value from MR05). When we compare 
our results with the only other published system parameters by 
109, we find that 109 obtain from their combined photometric and 
spectroscopic analysis a lower value of the eccentricity (e=0.19). 
When we fit their data using our RV solution method, we obtain 
e = 0.25, O ^ 47° and To = -0.4. However, when we use 
Q. = 80° and To = -0.1 as fixed values we obtain a^- 2 value that 
increases by less than 1 (i.e. less than l<x), indicating that this 
solution is also plausible with their data. The average y = -28.0 
km s _1 is consistent with the cluster radial velocity of -25.0 km 
s _1 (Valitova et al. 1990): providing extra evidence of the sys- 
tem's membership to the cluster. 

The derived value of i = 73° indicates that if we aim to de- 
termine precise masses, we must accurately determine the values 
of Ka and Kb- Given that we obtain a range of about 6 km s _1 
in both Ka and Kb this implies that the minimum mass error is 
already on the order of a few tenths (if we just use the analytical 
calculation and the inferred inclination). The resulting values for 
Msin 3 i are listed in Table [2] where it can be seen that the vari- 
ation is of the order of 0.3 M (all however within 3cr of each 
other). In Sect. [4] we already discussed the sensitivity of each of 
the methods in obtaining the RVs, and we would like to point 



4 MOOG was developed by Chris Sneden, see: 
http ://verdi . as . utexas . edu/moog .html 
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out the difficulty of attempting to model the type of stellar sys- 
tem considered here. What hampers the precise determination of 
the RVs of each component is the significant broadening of the 
spectral line profiles due to the rotation of both stellar compo- 
nents, and the superimposed pulsations. 

The largest differences in the results come from using the 
same method (IRAF) but choosing a different cut-off region in 
the CCF profiles to model. While the two-Gaussian function 
did not show reasonable results for IRAFfull, an adequate 3- 
Gaussian fitting to the CCF seemed to fit the profile better in- 
deed. This third low amplitude component could be a sign of a 
third body in the orbit. 

When combining the photometric and spectroscopic data 
(Sectj5]l it was difficult to arrive at a solution that (iteratively) 
converged using the parameters listed in Table [3] as free param- 
eters. This is most likely a result of the rapid rotation causing 
surface temperature variations at different latitudes (sometimes 
a difference of up to 1,500 K for very rapid rotation). The photo- 
metric light curve is distorted due to these variations, and these 
effects are not taken into account in the modeling (the assump- 
tion with PHOEBE is that the star rotates as a rigid body, without 
differential rotatior0). Additionally, one of the components has 
visible photometric variability on time scales of ~1 hour (MR05, 
Amado et al. 2006, Costa et al. 2007), and although the light 
curves were binned to reduce the effects of pulsation, there are 
nonetheless effects that can not be eliminated, such as those dur- 
ing eclipse. 

We obtain component masses of 1 .78 + 0.24 and 1 .70 + 0.22 
M for primary and secondary stars respectively, values that are 
lower by ~ 0.2 M than those from 109 (2.06 and 1 .87 M Q ). Our 
results are consistent with the values obtained by Mermillod & 
Mayor (1990) who, through a study of red giants, determined the 
main sequence (MS) turn-off mass of IC 4756 to be between 1.8 
and 1.9 M Q . This would also be consistent with the hypothesis of 
the more massive component beginning to turn off the MS, while 
also implying that the secondary star is slightly hotter than the 
primary. 

We obtain primary and secondary radius of 4.03 +0.11 and 
2.37 + 0.07 R Q , respectively, luminosities of 52.42 ± 2.85 and 
22.23 ± 1.33 L e , and logg values of 3.48+0.08 and 3.92+0.08 
(see Table [3). These values are lower than those quoted by 109, 
due primarily to the lower fitted masses we obtain. 

We finally fixed the values of r eff at 7500-7750 and 8000- 
8250 K respectively. This choice is justified by the following 
reasons: (a) at an earlier stage during this work, modeling the 
photometric light curve produced similar answers (7,560 and 
8,030 K) while fixing various other parameters, (b) modeling the 
light curve while fixing the T e g at other values results in worse 
fits or non-convergence, (c) the photometric Stromgren indices 
are consistent with these ranges of values (Sect. |5), and (d) the 
spectroscopic determination of T e ff from the LSD profiles and 
synthetic spectra also arrives at these values (Sect|7J. 

The primary (hotter and more massive) star is the pulsating 
component that clearly shows line-profile variations. Our data 
do not allow us to determine if the secondary component is also 
pulsating as the SNR is too low to detect low-amplitude variabil- 
ity, if any. The results in Sect. [5] and Table [3] do not discard the 
possibility of two pulsating components — both stars lie within 
the 5 Scuti instability strip close to the blue edge (see Fig. 13 in 
Pamyatnykh 2000). 
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Fig. 7. Abundances of the primary (filled circles) and secondary 
(open circles) stars compared to the adopted solar values. The 
bars represent errors of ± 1 cr. 



We have explored various forms of determining the RVs 
from the 70 spectra collected between 2005 and 2007 as well as 
using photometric observations spanning over 3000 days. Taking 
into account the uncertainties, our results are consistent with the 
109 ones, however, apart from the RV dataset derived with the 
IRAF method, we obtained lower masses and this would imply 
a downward revision of these values to 1.8 and 1.7 M . 



9.2. Abundances 

The resulting abundances given in Table H indicate for the pri- 
mary component a sub-solar metal value of [Fe/H] = -0.28 
consistent with the values of -0.15 and -0.22 determined for 
the cluster IC 4756 derived from samples of F-G single star 
cluster members by Jacobson et al. (2007) and Thogersen 
et al. (1993), respectively. The results of 109 also showed 
that this system component's global parameters fitted better 
to a sub-solar metallicity evolutionary track (Z ~ 0.008, 
where Z is the initial metal mass fraction). For HD 172189 
we obtain the following abundances (primary/secondary): 
[Na/Fe] = +0.03/-, [Al/Fe] = +0.00/-, [Si/Fe] = +0.08/+0.00, 
[Ca/Fe] = +0.03/+0.04, and [Ni/Fe] = +0.02/+0.04. These val- 
ues are consistent with those of IC 4756 according to Jacobson 
et al. (2007). However, their values are higher than ours by up to 
-0.3 dex (they quote +0.6, +0.3, +0.34, +0.07, +0.08, respec- 
tively). This could be in part due to the adopted log —gf values 
(Jacobson et al. 2007). 

The abundances relative to hydrogen for the primary star 
are: [Fe/H] = -0.28, [Si/H] = +0.3, [Ca/H] = -0.03, 
[Ni/H] = -0.10, [Na/H] = +0.04. These values are compara- 
ble to Jacobson et al. (2007), who obtained -0.15, +0.19, -0.08, 
-0.07, +0.^3 for IC 4756. The results for the secondary star are 
[Fe/H] = +0.4, [Si/H] = +0.4, [Ca/H] = +0.6, [Ni/H] = +0.6, on 
average ~0.5 dex higher than the primary star, resulting in en- 
hanced abundances compared to the Sun. Fig.|7]shows the abun- 
dances relative to the solar ones. The enhanced abundances of 
the secondary star could also be an effect of modeling a spec- 
trum with an inadequate normalisation from the disentangling. 
Higher SNR disentangled spectra resulting from the full orbital 
phase coverage would be able to confirm these different abun- 
dances. 



5 See the P HOEBE manual, 
http://phoebe.fiz. uni-lj .si/ 



available 



at 



Here, Jacobson et al. used the gf values from Luck (1994). 
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10. Conclusions 

We analysed a set of spectroscopic data obtained between 2005 
and 2007 with the aim of determining the absolute compo- 
nent masses, radii, luminosities, and effective temperatures of 
the system, as well as undertaking the first abundance analy- 
sis of the HD 172189 disentangled spectra. From a combined 
photometric and radial velocity analysis, we derive for the pri- 
mary and secondary stars, respectively, masses of 1.78+0.24 and 
1.70+0.22 M radii of 4.03+0.11 and 2.37+0.07 R , luminosi- 
ties of 52.42+2.85 and 22.23+1.33 L , and logg of 3.48+0.08 
and 3.92+0.08 dex. The analysis of the spectra indicates a T e g 
of 7600+150 K and 8100+150 K. From the RV data, we derive 
the system orbital parameters: e = 0.29 + 0.02, Q = 77° ± 5 
(from KOREL), and y = -28.02 ± 0.37 km s 1 (from other RV 
methods), and from the photometric data i = 73.2° + 0.6 and 
a = 20.32 ± 0.06 R Q - The spectroscopic data also confirm the 
orbital period of the system: 5.70198 days, as derived by MR05. 
We also obtain a spectral type of A6V-A7V based on Stromgren 
photometry. 

We applied four methods to calculate the RVs and obtained 
different results for each. We subsequently showed the limita- 
tions of each technique for analysing a system with two rapidly- 
rotating components, one of these having the additional com- 
plications of showing significant profile-variations due to pul- 
sations. Obtaining higher SNR spectra with continuous phase 
coverage would help to determine the systematic differences be- 
tween each of the methods. 

We have disentangled the spectra of both components and 
determined the rotational velocities: vskn = 78+3 km s _1 and 
74±4 km s _1 . Coupling these values with the results for R and 
; we obtain rotational periods P rot ,A = 2.50 ± 0.10 days and 
^rot.B = 1.55 ± 0.09 days. We subsequently derived a metallicity 
of [Fe/H] = -0.28 dex and abundances of [Si/H] = +0.3, [Ca/H] 
= -0.03, [Ni/H] = -0.10 and [Na/H] = +0.04 for the primary 
star. These are consistent with the results published for the clus- 
ter IC 4756 by Jacobson et al. (2007). The sub-solar metallicity 
is also consistent with the findings from 109. 

The rapid rotation of both components, the non-synchronous 
rotation (based on the orbital period), the eccentricity, and the 
likely membership of IC 4756 suggests that the system is still 
detached. Based on the inferred age of the system of 0.8-0.9 
Gyr (Alcaino 1965, Mermillod & Mayor 1990) we estimate that 
the 1.8 M primary star is moving towards the end of its MS 
lifetime, when it begins to cool down and increase in luminos- 
ity. Mermillod & Mayor (1990) also determine that the main se- 
quence turn-off mass for this cluster is 1 .8-1 .9 M , which agrees 
with this hypothesis. 

In order to study the oscillations of the primary star, re- 
cently CoRot (Baglin et al 2006a,b, Michel et al. 2008) observed 
HD 172189 as a primary asteroseismic target. Our analysis of 
this data has given a thorough insight into the nature of this ob- 
ject, and hence forms a solid foundation for the subsequent as- 
teroseismic analysis. We look forward to the unravelling of the 
mysteries associated with this object, and to learning about the 
effects of rotation on the interior structure of the star from the 
pulsation frequencies. 
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